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1. Introduction 

A classification of random-matrix ensembles by sym- 
metries was first established by Dyson, in an influential 
1962 paper with the title "the threefold way: algebraic 
structure of symmetry groups and ensembles in quan- 
tum mechanics". Dyson's threefold way has since be- 
come fundamental to various areas of theoretical phy- 
sics, including the statistical theory of complex many- 
body systems, mesoscopic physics, disordered electron 
systems, and the field of quantum chaos. 

Over the last decade, a number of random-matrix 
ensembles beyond Dyson's classification have come to 
the fore in physics and mathematics. On the physics 
side these emerged from work on the low-energy Dirac 
spectrum of quantum chromodynamics, and from the 
mesoscopic physics of low-energy quasi-particles in dis- 
ordered superconductors. In the mathematical research 
area of number theory, the study of statistical correla- 
tions in the values of Riemann zeta and similar func- 
tions has prompted some of the same generalizations. 

In this article, Dyson's fundamental result will be re- 
viewed from a modern perspective, and the recent ex- 
tension of Dyson's threefold way will be motivated and 
described. In particular, it will be explained why sym- 
metry classes are associated with large families of sym- 
metric spaces. 



2. The framework 

Random matrices have their physical origin in the 
quantum world, more precisely in the statistical the- 
ory of strongly interacting many-body systems such as 
atomic nuclei. Although random-matrix theory is nowa- 
days understood to be of relevance to numerous areas of 
physics - see e.g. Guhr's article in this volume - quan- 
tum mechanics is still where many of its applications 



lie. Quantum mechanics also provides a natural frame- 
work in which to classify random-matrix ensembles. 

Following Dyson, the mathematical setting for clas- 
sification consists of two pieces of data: 

• A finite-dimensional complex vector space "V with 
a Hermitian scalar product (•,•), called a unitary 
structure for short. (In physics applications, "V 
will usually be the truncated Hilbert space of a 
family of quantum Hamiltonian systems.) 

• On y there acts a group G of unitary and anti- 
unitary operators (the joint symmetry group of the 
multi-parameter family of quantum systems). 

Given this setup, one is interested in the linear space 
of self-adjoint operators on "V - the Hamiltonians H - 
with the property that they commute with G. Such a 
space is reducible in general, i.e. the matrix of H de- 
composes into blocks. The goal of classification is to 
enumerate the irreducible blocks that occur. 

2.1. Symmetry groups. — Basic to classification is 
the notion of a symmetry group in quantum Hamilto- 
nian systems, a notion that will now be explained. 

In classical mechanics the symmetry group Go of a 
Hamiltonian system is understood to be the group of 
canonical transformations that commute with the phase 
flow of the system. An important example is the rota- 
tion group for systems in a central field. 

In passing from classical to quantum mechanics, one 
replaces the classical phase space by a quantum me- 
chanical Hilbert space Y and assigns to the symme- 
try group Go a (projective) representation by unitary C- 
linear operators on "V . Beside the one-parameter con- 
tinuous subgroups, whose significance is highlighted 
by Noether's theorem, the components of Go not con- 
nected with the identity play an important role. A promi- 
nent example is provided by the operator for space re- 
flection. Its eigenspaces are the subspaces of states with 
positive and negative parity, which reduce the matrix of 
any reflection-invariant Hamiltonian to two blocks. 

Not all symmetries of a quantum mechanical system 
are of the canonical, unitary kind: the prime counterex- 
ample is the operation of inverting the time direction, 
called time reversal for short. In classical mechanics 
this operation reverses the sign of the symplectic struc- 
ture of phase space; in quantum mechanics its algebraic 
properties reflect the fact that inverting the time direc- 
tion, t \— > —t, amounts to sending i = \f^-Y to — i. In- 
deed, time t enters in the Dirac, Pauli, or Schrodinger 
equation as ihd/dt. Therefore, time reversal is repre- 
sented in the quantum theory by an anft'-unitary opera- 
tor T, which is to say that T is complex anti-linear: 

T(zy)=zTy (zeC, yeV), 
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and preserves the Hermitian scalar product or unitary 
structure up to complex conjugation: 

(Txj/ u Txi/ 2 ) = (YuYi) = (V2,Vi) • 

Another operation of this kind is charge conjugation in 
relativistic theories such as the Dirac equation. 

By the symmetry group G of a quantum mechanical 
system with Hamiltonian H, one then means the group 
of all unitary and anti-unitary transformations g of Y 
that leave the Hamiltonian invariant: gHg~ l = H. We 
denote the unitary subgroup of G by Go, and the set 
of anti-unitary operators in G by G\ (not a group). If 

Y carries extra structure, as will be the case for some 
extensions of Dyson's basic scheme, the action of G on 

Y has to be compatible with that structure. 

The set G\ may be empty. When it is not, the com- 
position of any two elements of G\ is unitary, so every 
g € G\ can be obtained from a fixed element of G\, say 
T, by right multiplication with some U G Go: g = TU. 
In other words, when G\ is non-empty the coset space 
G/Go consists of exactly two elements, Go and T Gq = 
G\. We shall assume that T represents some inversion 
symmetry such as time reversal or charge conjugation. 
T must then be a (projective) involution, i.e. T 2 = z, x Id 
with z a complex number of unit modulus, so that con- 
jugation by T 2 is the identity operation. Since T is 
complex anti-linear, the associative law T 2 T = T T 2 
forces z to be real, and hence T 2 = ±Id. 

Finding the total symmetry group of a Hamiltonian 
system need not always be straightforward, but this com- 
plication will not be an issue here: we take the symme- 
try group G and its action on the Hilbert space Y as 
fundamental and given, and then ask what is the corre- 
sponding symmetry class, meaning the linear space of 
Hamiltonians on Y that commute with G. 

For technical reasons, we assume the group Go to be 
compact; this is an assumption that covers most (if not 
all) of the cases of interest in physics. The non-compact 
group of space translations can be incorporated, if nec- 
essary, by wrapping the system around a torus, whereby 
translations are turned into compact torus rotations. 

While the primary objects to classify are the spaces 
of Hamiltonians H, we shall focus for convenience on 
the spaces of time evolutions U t = e~ ltH / n instead. This 
change of focus results in no loss, as the Hamiltonians 
can always be retrieved by linearizing in t at t = 0. 

2.2. Symmetric spaces. — We appropriate a few ba- 
sic facts from the theory of symmetric spaces. 

Let M be a connected m-dimensional Riemannian 
manifold and p a point of M. In some open subset N p 
of a neighborhood of p there exists a map s p :N p ^N p , 
the geodesic inversion with respect to p, which sends 



a point x 6 N p with normal coordinates (x\,...,x m ) to 
the point with normal coordinates (— x\ —x m ). The 
Riemannian manifold M is called locally symmetric if 
the geodesic inversion is an isometry, and is called glob- 
ally symmetric if s p extends to an isometry s p :M^M, 
for all p G M. A globally symmetric Riemannian man- 
ifold is called a symmetric space for short. 

The Riemann curvature tensor of a symmetric space 
is covariantly constant, which leads one to distinguish 
between three cases: the scalar curvature can be posi- 
tive, zero, or negative, and the symmetric space is said 
to be of compact type, Euclidean type, or non-compact 
type, respectively. (In mesoscopic physics each type 
plays a role: the first provides us with the scattering ma- 
trices and time evolutions, the second with the Hamil- 
tonians, and the third with the transfer matrices.) The 
focus in the current article will be on compact type, as 
it is this type that houses the unitary time evolution op- 
erators of quantum mechanics. The compact symmetric 
spaces are subdivided into two major subtypes, both of 
which occur naturally in the present context, as follows. 

2.3. Type II. — Consider first the case where the anti- 
unitary component G\ of the symmetry group is empty, 
so the data are (Y,G) with G = G . Let fy(Y) de- 
note the group of all complex linear transformations 
that leave the structure of the vector space Y invariant. 
Thus ffl (Y) is a group of unitary transformations if Y 
carries no more than the usual Hermitian scalar prod- 
uct; and is some subgroup of the unitary group if Y 
does have extra structure (as is the case for the Nambu 
space of quasi-particle excitations in a superconductor). 
The symmetry group Go, by acting on Y and preserving 
its structure, is contained as a subgroup in & (Y). 

Let now H be any Hamiltonian with the prescribed 
symmetries. Then the time evolution 1 1— > U t = e~ ltH / h 
generated by H is a one-parameter subgroup of 'W (Y) 
which commutes with the Go-action. The total set of 
transformations U t that arise in this way is called the 
(connected part of the) centralizer of Go in <%f(Y), and 
is denoted by Z. This is the "good" set of unitary time 
evolutions - the set compatible with the given symme- 
tries of an ensemble of quantum systems. 

The centralizer Z is obviously a group: if U and V 
belong to Z, then so do their inverses and their product. 
What can one say about the structure of the group Z? In 
essence, this question was answered by H. Weyl in his 
famous treatise "The Classical Groups". 

Since Go is compact by assumption, its group action 
on Y is completely reducible, and Y is guaranteed to 
have an orthogonal vector space decomposition 

X X 
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where the sum is over (classes of equivalent) irreducible 
Go-representations A , the Vx are irreducible representa- 
tion spaces for the Go-action, and mx is the multiplic- 
ity of occurrence of the representation type A. (Go acts 
trivially on C m K) The subspaces y x ~ V x <8> C m >- will be 
called the Go-isotypic components of y. For example, 
if Go is the rotation group SO(3), the Go-isotypic com- 
ponents of y are the subspaces of states with definite 
total angular momentum, say A =0,1,2,...; Vx then 
is an SO (3) -irreducible representation space of dimen- 
sion 2A + 1 ; and mx is the number of times a multiplet 
of states with angular momentum A occurs in "V . 

Now consider any U G Z. Since U commutes with 
the Go-action, it does not connect different Go-isotypic 
components. (Indeed, in the example of SO(3)-invariant 
dynamics, angular momentum is conserved and transi- 
tions between different angular momentum sectors are 
impossible.) Thus every Go-isotypic component fx is 
an invariant subspace for the action of Z on "V , and Z 
decomposes as Z = Ylx %x with blocks Zx = Z| y . But 
one can say even more: because U G Z commutes with 
Go and Vx is Go-irreducible, U must act like the iden- 
tity on Vx by Schur's lemma. Therefore, Zx acts non- 
trivially only on the factor C'" 1 in Y x - Vx ® C mi , so 

Z = \\Zx — \\Vi(mx) (direct product) , 

X X 

if "V is a unitary vector space with no extra structure. In 
the presence of extra structure (which, by compatibility 
with the Go-action, restricts to every subspace the 
factor Zx is some subgroup of U(m^). In all cases, Z is 
a direct product of connected compact Lie groups Zx . 

The focus now shifts from Z to any one of the Zx- 
So we fix some M := Z^. Since M is a group, the op- 
eration of taking the inverse, U > U~ l , makes sense 
for all U G M. Moreover, being a compact Lie group, 
the manifold M admits a left- and right-invariant Rie- 
mannian structure in which the inversion U i-> U~ l is 
an isometry. By translation one gets an isometry 

for every U\ G M. All these maps su l are globally de- 
fined, and the restriction of sy l to some neighborhood of 
U\ coincides with the geodesic inversion w.r.t. U\. Thus 
M is a symmetric space by the definition given above. 
Symmetric spaces of this kind (with group structure) 
are called type II. 

2.4. Type I. — Consider next the case G\ / 0, where 
some anti-unitary symmetry T is present. As before, let 
Z be the connected component of the centralizer of Go 
in W(y). Conjugation by T, 



is an automorphism of (Y) and, owing to T 2 = ±Id, 
T is involutive. Because Go C G is a normal subgroup, 
T restricts to an involutive automorphism (still denoted 
by t) of Z. Now recall that T is complex anti-linear 
and the good Hamiltonians are subject to THT 1 = H. 
The good time evolutions U t = e~ ltH / n clearly satisfy 
z(U t ) = U- t = Uf~ l . Thus the good set to consider is 

Jl :={U eZ\U = r(uy 1 } . 

is a manifold, but in general is not a Lie group. 

Further details depend on what t does with the fac- 
torization Z = Ylx Zx- If V\ is a Go-isotypic component 
of "f, then so is TYx, since T normalizes Go- Thus 
TYx = y~x for some representation type A. If A / A, 
the involutive automorphism T just relates U G Z^ with 
t(U) GZr , whence no intrinsic constraint on Zx results, 
and the time evolutions (U,T(U)~ r ) EZxxZ-^ consti- 
tute a type-II symmetric space, as before. 

A novel situation occurs when A = A, in which case 
T maps the group Zx onto itself. Let therefore A = A , 
put K = Zx for short, and consider 

M :={U £K\U = r(uy 1 } . 

Note that if two elements p,po of K are in M, then so is 
the product poP~ l po- The group K acts on M C K by 

k-U = kUx{ky l (k€K), 

and this group action is transitive, i.e. every U G M can 
be written as U = kz(k)~ l with some k G K. (Finding k 
for a given U is like taking a square root, which is pos- 
sible since exp : LieK —> K is surjective.) There exists 
a i^-invariant Riemannian structure for M such that for 
all po G M the mapping s Po :M^M defined by 

s po(p) =P0P~ l P0 , 

is the geodesic inversion w.r.t. po G M. Thus in this nat- 
ural geometry M is a globally symmetric Riemannian 
manifold and hence a symmetric space. The present 
kind of symmetric space is called type I. If K x is the 
set of fixed points of % in K, the symmetric space M is 
analytically diffeomorphic to the coset space K/K x by 

K/K t ^M CK , UK X ^ Uz{Uy l , 

which is called the Carton embedding of K/K r into K. 

In summary, the solution to the problem of finding 
the unitary time evolution operators that are compatible 
with a given symmetry group and structure of Hilbert 
space, is always a symmetric space. This is a valuable 
insight, as symmetric spaces are rigid objects and have 
been completely classified by Cartan. 

If we keep the dimension of y variable, the symmet- 
ric spaces that occur must be those of a large family. 
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3. Dyson's threefold way 

Recall the goal: given a Hilbert space Y and a sym- 
metry group G acting on it, one wants to classify the 
(irreducible) spaces of time evolution operators U that 
are "compatible" with G, meaning 

V =g Ug 1 =g l U- 1 g^ 1 (for all*, G G a ) ■ 

As we have seen, the spaces that arise in this way are 
symmetric spaces of type I or II depending on the nature 
of the time reversal (or other anti-unitary symmetry) T. 

An even stronger statement can be made when more 
information about the Hilbert space Y is specified. In 
Dyson's classification, the Hermitian scalar product of 
Y is assumed to be the only invariant structure that ex- 
ists on Y. With that assumption, only three large fami- 
lies of symmetric spaces arise; these correspond to what 
we call the Wigner-Dyson symmetry classes. 

3.1. Class A. — Recall that in Dyson's case, the con- 
nected part of the centralizer of Go in W{Y) is a direct 
product of unitary groups, each factor being associated 
with one Go-isotypic component Yx of Y. The type-II 
situation occurs when the set G\ of anti-unitary sym- 
metries is either empty or else exchanges different Y\ . 
In both cases, the set of good time evolution operators 
restricted to one Go-isotypic component Yx is a unitary 
group U(mi), with mx being the multiplicity of the ir- 
reducible Go-representation X in Yx . 

The unitary groups U(N = m^ ) or to be precise, their 
simple parts SU(Af), are called type-II symmetric spaces 
of the A family or A series - hence the name class A. 
The Hamiltonians H, the generators of time evolutions 
U t = &~ ltH l Tl , in this class are represented by complex 
Hermitian N x N matrices. By putting a U(A r )-invariant 
Gaussian probability measure 

exp (-Tr// 2 /2a 2 ) dH (a G R) 

on that space, one gets what is called the GUE - the 
Gaussian Unitary Ensemble - which defines the Wigner- 
Dyson universality class of unitary symmetry. 

3.2. Classes AI and All. — Consider next the case 
G\ 0, with anti-unitary generator T. Let Yx be any 
Go-isotypic component which is invariant under T (the 
type-I situation). The mapping U \— > TUT~ l = t{U) 
then is an automorphism of the groups \J(Yx), Go and 
K = Zx— U(m^). If K t is the subgroup of fixed points 
of T in K, the space of good time evolutions can be 
identified with the symmetric space K/K x by the Cartan 
embedding. The task is to determine K t . As was em- 
phasized by Dyson, the answer for K t does not follow 
from any single piece of data, but is determined by the 
combination of three anti-unitary involutions on Yx . 



The first of these is the standard operation of taking 
the complex conjugate (w.r.t. the complex structure of 
Yx), denoted by y h-> y as usual. Recall that the uni- 
tary vector space Yx decomposes as an orthogonal sum 
of mx identical copies of an irreducible representation 
space Vx for the compact group Go: 

(1) Yx = ©Vf©...® V^ -V x ® <£ mx ■ 

From the multitude of such decompositions, we select 
one that is invariant under complex conjugation: if y 
lies in a subspace Vv , so does y. 

Next we look at the operation go h-> g on any one 
of the Go-irreducible subspaces of Yx, say Vx = vj. . 
Since X = X by assumption, the Go-action on Vx is uni- 
tarily equivalent to its complex conjugate. Thus there 
exists some unitary transformation s G U(V^) such that 

go = s~ 1 gos 

holds for every go G Go. Given s, one defines an anti- 
unitary operator S on Yx by 

S(v^ + ...+v {m ^) =^W + ...+j^KT, 

where the decomposition Q is invoked. By construc- 
tion, S commutes with the action of Go- Since the Go- 
action on Vx is irreducible, Schur's lemma applied to S 2 
forces S to be a projective involution, and associativity 
(S 2 ■ S = S ■ S 2 ) results in S 2 = ±\ n . 

Recall next that go | — ► TgoT^ 1 is an automorphism 
of Go, and remains so when the Go-action is restricted 
to Yx- Since all of the Go-representations in the iso- 
typic component Yx are equivalent, there exists a uni- 
tary transformation R G V>{Yx) such that 

Tg Q T- 1 =R- l g R (for all go G Go) 

holds as an operator identity on Yx ■ Note that the com- 
position RT intertwines Go-actions: RTgo = goRT, but 
changes the complex structure of Yx (by anti-linearity 
of T). A better object to consider is RT o 5 which, be- 
ing composed of two anti-unitary operators, is unitary. 
RTS commutes with the Go-action and thus lies in the 
centralizer K. Using it, one defines another anti-unitary 
operator T on Yx and an automorphism z' of \J(Yx) by 

T'y = RTSy , z'(k) = T'kT'~ l . 

T determines a complex bilinear form Q on Yx by 

60i . ¥2) = {T'W\ , W2) n (for all fi.^e^). 

The remaining steps toward identifying K t depend 
on the nature of R. Consider first the easy case where 
the automorphism X of Go is inner, i.e. R G Go- Then 
k = %(k) for k G K is equivalent to k = RTk(RT)-\ 
which in turn amounts to k = RTSk(RTS)^ i = t'(k). 
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Hence another description of K z is to say that its ele- 
ments k are the unitary transformations of Yi that cen- 
tralize Go and leave the pairing Q invariant (the latter 
follows from T'k = IcT' and invariance of (•,•)■%). 

By iterating k = T'(k) = T ,2 (k) one infers that T' 2 
commutes with the ^-action on "V%. But T' 2 also com- 
mutes with T' and with the Go-action, which implies 
T' 2 = e x ly x with £ = ±1, by standard reasoning. From 

Q(Vi,V2) = {T' 2 y x ,T'Vi)v x =eQ{V2,Vi) 

one sees that the pairing Q is symmetric for e = +1, 
and skew for e = — 1. The Lie group K z is now easily 
identified. Since the unitary operator RTS is an element 
of K ~ U(m^), Q restricts to a pairing on the factor C mi 
in the decomposition Yx — Vx ®C m;l . Thus for e = +1, 
K x can be viewed as a subgroup of U(m^) that preserves 
a symmetric pairing (or orthogonal structure) on C mx ; 
consequently K r ~ O(m^). For e = —1, the multiplicity 
mi must be even, and K x preserves a skew pairing (or 
symplectic structure); in that case K x ~ USp(m^), the 
unitary symplectic group. 

In the general case (R ^ Go) drawing these conclu- 
sions is more difficult, and one must exploit the rigidity 
of the symmetric space K/K r under deformations of the 
automorphism t(U) = TUT^ 1 . Actually, the method 
used by Dyson to handle the general case is quite dif- 
ferent: Dyson chooses to regard as a real vector 
space (with complex structure) and expresses the group 
action of G by orthogonal matrices. In this real setup 
he then exploits a deep theorem of Weyl on the struc- 
ture of group algebras and their commutator algebras, 
which leads him to the conclusion that the above two 
possibilities are in fact the only ones that can occur in 
the present type-I situation. Thus there is a dichotomy 
for the sets of good time evolutions M ~K/K Z : 

Class AI : K/K z ~ U(N)/0(N) (N = m k ), 

Class All : K/K T ~ \]{2N) /USp(2A0 (2N = m x ) . 

Again we are referring to symmetric spaces by the names 
they - or rather their simple parts S\J{N)/SO(N) and 
SU(2A0/USp(2A0 - have in the Cartan classification. 
Be warned that Weyl's theorem by itself does not allow 
to decide between the alternatives AI or All for a given 
data set (Y x ,Go,T) (rather, to do so you must determine 
the "Wigner type" of the G-representation on Yx). 

In the case where the Go-automorphism T is inner 
(which actually covers most of the known examples of 
physical interest) Dyson's reasoning is basically identi- 
cal to the one reviewed above. There, as we have seen, 
the dichotomy is ruled by the number e computed from 
(J 1 ) 2 = RTS RTS = £xlf r Important examples are 



provided by physical systems with spin-rotation sym- 
metry, Go = SU(2), and time-reversal symmetry. The 
physical operation of time reversal, T, commutes with 
spin rotations, so T is inner here with R = 1 . On states 
with spin one has T 2 = (-l) 2 ^ and S 2 = (-l) 2|s| , 
which gives T' 2 = +1 in all cases. Thus time-reversal 
invariant systems with no symmetries other than energy 
and spin are always class AI. By breaking spin-rotation 
symmetry (Go = {Id}, so S 2 = 1) while maintaining T- 
symmetry for states with half-integer spin (say single 
electrons, which carry spin \S\ = 1/2), one gets T' 2 = 
T 2 = — 1 , thereby realizing class All. 

The Hamiltonians H, obtained by passing to the tan- 
gent space of K/K r at unity, are represented by Hermi- 
tian matrices with entries that are real numbers (class 
AI) or real quaternions (class All). If you put ^-invari- 
ant Gaussian probability measures on these spaces, you 
get the Wigner-Dyson universality classes of orthogo- 
nal resp. symplectic symmetry. In mesoscopic physics 
these are realized in disordered metals with time-reversal 
invariance (absence of magnetic fields and magnetic im- 
purities). Spin-rotation symmetry is broken by strong 
spin-orbit scatterers such as gold impurities. 



4. Disordered superconductors 

When Dirac first wrote down his famous equation 
in 1928, he assumed that he was writing an equation 
for the wavefunction of the electron. Later, because of 
the instability caused by negative-energy solutions, the 
Dirac equation was reinterpreted (via second quantiza- 
tion) as an equation for the fermionic field operators of 
a quantum field theory. A similar change of viewpoint 
is carried out in reverse in the Hartree-Fock-Bogoliubov 
mean field description of quasi-particle excitations in 
superconductors. There, one starts from the equations 
of motion for linear superpositions of the electron cre- 
ation and annihilation operators, and reinterprets them 
as a unitary quantum dynamics for what might be called 
the quasi-particle "wavefunction". 

In both cases - the Dirac equation and the quasi- 
particle dynamics of a superconductor - there enters a 
structure not present in the standard quantum mechan- 
ics underlying Dyson's classification: the field opera- 
tors for fermionic particles are subject to a set of re- 
quirements called the canonical anti-commutation rela- 
tions, and these are preserved by the quantum dynam- 
ics. Therefore, whenever second quantization is undone 
(assuming it can be undone) to return from field opera- 
tors to wavefunctions, the wavefunction dynamics is re- 
quired to preserve some extra structure. This puts a lin- 
ear constraint on the allowed Hamiltonians H. For our 
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purposes, the best viewpoint to take is to attribute the 
extra invariant structure to the Hilbert space i / , thereby 
turning it into a Nambu space. 

4.1. Nambu space. — Starting from the standard for- 
malism of second quantization, consider a set of single- 
particle creation and annihilation operators c« and c a , 
where a = l,...,N labels single-particle states that are 
orthogonal to each other. Such operators are subject to 
the canonical anti-commutation relations 

(2) c lt c p + c fi c a = <5aj3 > 

C «4 + C j3 C " = = CaC P + C P Ca ■ 

When written in terms of c a + c« and i(c« — c a ), these 
become the defining relations of a Clifford algebra. Field 
operators are linear combinations y = £ a {u a c* a + v a c a ) 
with complex coefficients u a and v a . 

Now take H to be some Hamiltonian which is quadratic 
in the creation and annihilation operators: 

H = Y* A aP C a C P + ( B «/? c «4 + S ap c fl c a) , 
aft a/3 

and let H act on field operators y by the commutator: 
H • y = [H,y]. The time evolution of y is then deter- 
mined by the equation 



which integrates to y{t) = er itH l h ■ y(0), and is easily 
verified to preserve the relations ©. 

The dynamical equation © is equivalent to a sys- 
tem of linear differential equations for the amplitudes 
u a and v a . If these are assembled into vectors, and the 
A a p and B a p into matrices, equation © becomes 

dt [yj " h \-B -A) \y) ■ 

The Hamiltonian matrix on the right-hand side has some 
special properties due to B a p = —Bp a (from c a cp = 
—cpc a ) and A a p = Ap a (from H being self-adjoint as 
an operator in Fock space). To keep track of these prop- 
erties while imposing some unitary and anti-unitary sym- 
metries, it is best to put everything in invariant form. 

Let V be the complex vector space of annihilation 
operators v = Y, a v a Ca, and view the creation operators 
u = H a u a cli as lying in the dual vector space V*. The 
field operators y = u + v then are elements of the direct 
sum V* © V =: "V , called Nambu space. On "f there 
exists a canonical unitary structure expressed by 

(mi + Vi , u 2 + v 2 ) = ^ ("la"2a + vi a v 2a ) . 
a 



A second canonical structure on 'f is given by the sym- 
metric C-bilinear form 

{ui + Vi , U 2 + V 2 } = Y ("l«V 2 a + U 2a Via) ■ 
a 

Note that {yi, y 2 } agrees with the anti-commutator of 
the field operators, y\ y 2 + y 2 yi, by the relations ©. 

Now recall that the quantum dynamics is determined 
by a Hamiltonian H that acts on y by the commutator 
H y = [H,y\. The one-parameter groups 1 1— ► e^ ltH / h 
generated by this action (the time evolutions) preserve 
the symmetric pairing: 

{^ 1 ,^ 2 } = {e-^.^,e- i ^-^ 2 }, 

since the anti-commutation relations © do not change 
with time. They also preserve the unitary structure, 

{yuy 2 ) = ^- mlh -w^- mlK W2), 

because probability in Nambu space is conserved. (Phys- 
ically speaking, this holds true as long as H is quadratic, 
i.e. many-body interactions are negligible.) 

One can now pose Dyson's question again: given 
Nambu space "V and a symmetry group G acting on it, 
what is the set of time evolution operators that preserve 
the structure of "V and are compatible with G? From 
Section |2] we know the answer to be some symmetric 
space, but which are the symmetric spaces that occur? 

4.2. Class D. — Consider a superconductor with no 
symmetries in its quasi-particle dynamics, so G = {Id}. 
(A concrete example would be a disordered spin-triplet 
superconductor in the vortex phase). The time evolu- 
tions U = e,~ ltH l n are then constrained only by invari- 
ance of the unitary structure and the symmetric pairing 
{-, ■} of Nambu space. These two structures are consis- 
tent; they are related by particle-hole conjugation C: 

{V1.V2} = (Cyi,y 2 ) , 

which is an anti-unitary operator with square C 2 = +Id. 
The condition {yi,y 2 } = {Uyi,Uy 2 } (invariance of 
an orthogonal structure) selects a complex orthogonal 
group, and imposing unitarity yields a real subgroup 
SO(X) ~ SO(4N) - a symmetric space of the D family. 

Since the time evolutions are a real orthogonal group, 
there exists a basis of "V (called Majorana fermions in 
physics) in which the matrix of \H G 50 {~f) is real skew, 
and that of H imaginary skew. The simplest random 
matrix model for class D, the SO-invariant Gaussian en- 
semble of imaginary skew matrices, is analyzed in the 
second edition of Mehta's book. From the expressions 
given by Mehta it is seen that the level correlation func- 
tions at high energy coincide with those of the Wigner- 
Dyson universality class of unitary symmetry. The level 
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correlations at low energy, however, show different be- 
havior defining a separate universality class. This uni- 
versal behavior at low energies has immediate physi- 
cal relevance, as it is precisely the low-energy quasi- 
particles that determine the thermal transport properties 
of the superconductor at low temperatures. 

4.3. Class Dili. — Let now magnetic fields and mag- 
netic impurities be absent, so that time reversal T is a 
symmetry of the quasi-particle system: G = {Id,T}. 
Following Section |2j the set of good time evolutions 
is M ~ K/K x with K = SO{f) and K T the set of fixed 
points of U t(Z7) = TUT- 1 in K. What is K z ? 

The time-reversal operator has square T 2 = —Id (for 
particles with spin 1/2), and commutes with particle- 
hole conjugation C, which makes Q := iCT a useful 
operator to consider. Since C by definition commutes 
with the action of K, and hence also with that of K % , the 
subgroup K x has an equivalent description as 

K T = {k e U(r) | k = QkQ- 1 = t(*)} . 

The operator Q is easily seen to have the following pro- 
perties: (i) Q is unitary, (ii) Q = Id, and (iii) Try Q = 
0. Consequently Q possesses two eigenspaces y± of 
equal dimension, and the condition k = QkQ~ l fixes a 
subgroup U(^+) x of \]{Y). Since Q contains a 

factor i = \J — 1 in its definition, it anti-commutes with 
the anti-linear operator T . Therefore the automorphism 
z exchanges U(^+) with U(^), and the fixed point set 
K T is the same as U(^+) ~ U(2N). Thus 

Af ~ K/K t ~ SO(4A0/U(2A0 , 

a symmetric space in the Dili family. Note that for par- 
ticles with spin 1/2 the dimension of "V + has to be even. 

By realizing the algebra of involutions C, T as C\ff = 
(ia*® l 2 /v)^and T\j/= (io y ® ^2n)V' tne Hamiltonians 
H in class Dili are brought into the standard form 

o). 

where the 2iV x 2Af matrix Z is complex and skew. 

4.4. Class C. — Next let the spin of the quasi-particles 
be conserved, as is the case for a spin-singlet super- 
conductor with no spin-orbit scatterers present, and let 
time-reversal invariance be broken by a magnetic field. 
The symmetry group of the quasi-particle system then 
is the spin-rotation group: G = Go = Spin(3) = SU(2). 

Nambu space "V can be arranged to be a tensor prod- 
uct V = W ® C 2 so that G acts trivially on W and by 
the spinor representation on the spinor space C 2 . Since 
two spinors combine to give a scalar, the latter comes 
with a skew-symmetric form e : C 2 x C 2 — > C. In a 
suitable basis, the anti-commutation relations © factor 



on particle-hole and spin indices. The symmetric bilin- 
ear form {•,•} of y correspondingly factors under the 
tensor product decomposition "V = W <g> C as 

{wi<g>Ji,w 2 <8>J2} = [wi,w 2 ] x e(s h s 2 ) , 

where [•, •] is a skew bilinear form on W, giving W the 
structure of a symplectic vector space. 

The good set M now consists of the time evolutions 
that, in addition to preserving the structure of Nambu 
space, commute with the spin-rotation group SU(2): 

M = {U € \3{y)\UC = CU,VR € SU(2) : RU = UR} . 

By the last condition, all time evolutions act trivially 
on the factor C 2 . The condition UC = CU, which ex- 
presses invariance of the orthogonal structure of Y, then 
implies that time evolutions preserve the symplectic pair- 
ing of W . Time evolutions therefore are unitary sym- 
plectic transformations of W , hence M = \JSp(W) ~ 
USp(2A^) - a symmetric space of the C family. The 
Hamiltonian matrices in class C have the standard form 

with A being Hermitian and B complex and symmetric. 

4.5. Class CI. — The next class is obtained by taking 
the time reversal T as well as the spin rotations 7? G 
SU(2) to be symmetries of the quasi-particle system. 

By arguments that should be familiar by now, the 
set of good time evolutions is a symmetric space M ~ 
K/K x with K = USp(#0 and K T the set of fixed points 
of T in K. Once again, the question to be answered is: 
what's K t ? The situation here is very similar to the one 
for class Dili, with W and USp(^) taking the roles 
of "¥ and SO(T). By adapting the previous argument 
to the present case, one shows that K t is the same as 
\]{W + ) ~ U(Af), where W + is the positive eigenspace 
of Q = iCT viewed as a unitary operator on W. Thus 

M~K/K^~ USp(2A0 /U(A0 . 

The standard form of the Hamiltonian matrices here is 

o) 

with the N xN matrix Z being complex and symmetric. 



5. Dirac fermions: the chiral classes 

Three large families of symmetric spaces remain to 
be implemented. Although these, too, occur in meso- 
scopic physics, their most natural realization is by 4d 
Dirac fermions in a random gauge field background. 
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Consider the Lagrangian L for the Euclidean space- 
time version of quantum chromodynamics with N c > 3 
colors of quarks coupled to an S\J(N C ) gauge field A^: 

L = iijf'f(d ll -A^yf+imyy . 

The massless Dirac operator D = ij^(d^ — An) anti- 
commutes with 75 = y y i y 2 y i . Therefore, in a basis 
of eigenstates of 75 the matrix of D takes the form 

If the gauge field carries topological charge v G Z, the 
Dirac operator D has at least |v| zero modes by the in- 
dex theorem. To make a simple model of the challeng- 
ing situation where A^ is distributed according to Yang- 
Mills measure, one takes the matrices Z to be com- 
plex rectangular, of size p x q with p — q = v, and puts 
a Gauss measure on that space. This random-matrix 
model for D captures the universal features of the QCD 
Dirac spectrum in the massless limit. 

The exponential of the truncated Dirac operator, q^ 
(where t is not the time), lies in a space equivalent to 
U(p + q) /U(p) x \J(q) - a symmetric space of the ALU 
family. We therefore say that the universal behavior of 
the QCD Dirac spectrum is that of symmetry class AIR. 

But hold on! Why are we entitled to speak of a sym- 
metry class here? By definition, symmetries always 
commute with the Hamiltonian, never do they anti-com- 
mute! (The relation D = —JsDy 5 is not a symmetry in 
the sense of Dyson, nor is it a symmetry in our sense.) 

5.1. Class AIII. — To incorporate the massless QCD 
Dirac operator into the present classification scheme, 
we adapt it to the Nambu space setting. This is done by 
reorganizing the 4-component Dirac spinor i/a, y as an 
8-component Majorana spinor *P, to write 

The 8x8 matrices are real symmetric besides sat- 
isfying the Clifford relations HT V +T v r^ = 2<5^ v . A 
possible tensor-product realization is 

l<g)C7 Z <g>l, r 1 = X ® Oy®Oy , 
T 2 = Oy ® Oy ® 1 , ^ = O z (g> Oy (g> Oy . 

The gauge field in this Majorana representation is stf^ = 

1 ® 1 ® (A { ~ ] -A^Oy) where Aj i ±) = \{A^ ±Aj l ) are 
the symmetric and skew parts of A^ £ su(N c ). 

The operator H = iP^d^, — =0^) is imaginary skew, 
therefore e ltH is real orthogonal. This means that there 
exists a Nambu space "f with unitary structure (-, •) and 
symmetric pairing {•,•}, both of which are preserved 
by the action of & ltH . No change of physical meaning or 
interpretation is implied by the identical rewriting from 



Dirac D to Majorana H. The fact that Dirac fermions 
are not truly Majorana is encoded in a U(l)-symmetry 
He i9Q = e, idQ H generated by Q = 1 ® 1 <g> a y . 

Now comes the essential point: since H obeys H = 
—H, the chiral "symmetry" H = —T=,HY=, with T=, = 
1 (g) a x Cg) 1 can be recast as a true symmetry : 

H = +r 5 HT 5 = THT 1 , 

with anti-linear T : *F 1— ► Y^. Thus the massless QCD 
Dirac operator is indeed associated with a symmetry 
class in the present, post-Dyson sense: that's class AIII, 
realized by self-adjoint operators on Nambu space with 
Dirac U(l)-symmetry and an anti-unitary symmetry T. 

5.2. Classes BDI and CII. — Consider Hamiltonians 
D still of the form © but now with matrix entries taken 
from either the real numbers or the real quaternions. 
Their one-parameter groups e ltD belong to two further 
families of symmetric spaces: 

Class BDI : SO(p + q) /SO{p) x SO(q) , 

Class CII : USp(2/> + 2q) /USp(2p) x USp(2^) . 

These large families are known to be realized as sym- 
metry classes by the massless Dirac operator with gauge 
group SU(2) (for BDI), or with fermions in the adjoint 
representation (for CII). For the details we must refer to 
Verbaarschot's paper, as there is no space left here. 
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